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Land Surface Interaction 

Robert E. Dickinson 

Land and Climate Modeling 

Historically, climate models have evolved out of general circula- 
tion models (GCMs) of the atmosphere, developed by meteorologists, 
and radiative convective models, developed by atmospheric physi- 
cists. These areas are my background. Land was thought of simply 
as a lower boundary condition to be put in as simply as possible. 
This was an implicit message in all of my academic training and was 
the basis for my first attempts to be involved in climate modeling. 
Conversely, many other disciplines have studied the role of soils and 
vegetation over land in great detail but traditionally with a very local 
focus, and with the atmosphere and radiation prescribed from 
observations rather than as part of the model. 

When GCMs were originally developed, requirements for accuracy 
and realism in the physical descriptions were much lower than they 
have become today. I recall conversations two decades ago where a 
model simulation of the tropical tropopause 40 degrees warmer than 
reality was regarded as satisfactory. Now only a few degrees is of 
considerable concern, a discrepancy small enough that it could 
result as readily from errors in the data themselves or their interpo- 
lation or the model layer structure as from errors in the model treat- 
ment of radiation. Except for the basic fluid equations, physical 
treatments in the original GCMs were all generally quite simple, rep- 
resented by a single equation and/or describable in one or two sen- 
tences. Land has been no exception. S. Manabe in 1969 first moved 
from the viewpoint of land as a boundary condition to land as an 
interactive part of the system with his “bucket" hydrology model; all 
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land points had a water-holding capacity with any surplus going 
into runoff, and otherwise the storage was determined by the bal- 
ance between precipitation or snowmelt and evapotranspiration. 
Surface albedos have been simply specified numerical values, often 
the same value for all land points, and the aerodynamic drag coeffi- 
cients {determined in principle by surface roughness) have had one 
value for land and one for ocean. 

Demands for better answers and increasing maturity of the sci- 
ence are driving climate models to much more detailed treatments of 
physical processes, much greater emphasis on validation of the 
assumed process physics, and in general greater complexity. These 
demands Eire shifting the requirements for training and the research 
emphases of individual scientists. Initial work with GCMs required a 
state-of-the-art knowledge of available methods for solving the 
model fluid equations and procedures for handling the large data 
sets generated by the models. Other than that, modelers needed 
only the ability to copy and implement the one-line descriptions of 
the various model physical processes. There was no urgent require- 
ment to validate or improve upon these descriptions, or to obtain 
observational data for model input or validation. Thus, it was possi- 
ble for a single highly skilled individual with the help of a couple of 
programmers to develop a GCM and carry out a research program 
with it. This traditional approach is still viable for programs empha- 
sizing improved understanding of atmospheric processes. However, 
the task of developing models with adequate realism and validation 
to meet the challenge of global change requires scientists to special- 
ize in particular aspects of a model, although still retaining at least 
some familiarity with the overall model structure and behavior. 
Indeed, we now need to develop teams or other collaborative group- 
ings of scientists to focus on the development of critical and still 
poorly developed aspects of climate models. Two particular examples 
that come to mind are the treatment of cloud-radiative interactions 
and land processes within models. 

Suppose you have decided to participate in a team effort to 
improve the treatments of land in climate system models with both 
the objective of relating this effort to development of this area as 
part of global change models and the immediate task of improving 
model projections of climate change. What must you know to get 
started? The issue is primarily the treatment of energy and water 
fluxes, since these have strong direct interactions with climate mod- 
els. At the same time, you must develop a framework for treatment 
of fluxes of carbon dioxide and other trace gases from the land sur- 
face. Currently, climate modelers start with scenarios of trace gas 
increases and do not attempt to make them interactive with the cli- 
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mate model. Practically, at present the future human contributions 
and natural feedbacks are too poorly known for there to be any 
practical improvement possible in model projections with inclusion 
of these feedbacks, but from a viewpoint of better understanding the 
system, they must be explored. These notes are an attempt to pro- 
vide the appropriate background material for someone who would 
like to start research in the treatment by climate models of surface 
water and energy fluxes. 

Our past knowledge base consists of two distinct threads of effort: 
(1) model sensitivity studies that give us better insight into the role 
of land processes and help sort out the relative importance of vari- 
ous factors, and (2) consideration of the detailed processes that 
must be modeled to represent land in climate models. To the extent 
that we are guided by the requirements of modeling the overall sys- 
tem, these two threads must proceed in parallel. In practice, we are 
working with a system of infinite complexity, and the level of 
abstraction vs. details required must be guided by overall modeling 
experience. For example, a decade ago we did not know enough 
about the treatment of land processes in climate models to be able 
to argue that we had to consider soils and vegetation as separate 
components. While there is still not general agreement as to this 
claim, a good case can be made that it is true. We can also show 
that the distinction between a forest and a grassland has a notice- 
able impact on a climate model, at least in the tropics. However, we 
are a long way from being able to distinguish through modeling the 
implications of switching between a maple and an oak forest, or pine 
and spruce, as might be provided by ecological modeling of response 
to climate change. Thus current research on feedbacks of surface 
fluxes with vegetation is emphasizing primarily short and intermedi- 
ate scales, that is, from model time steps of a few tens of minutes to 
the annual and interannual time scales. 


Sensitivity Studies 

A climate model with an atmospheric hydrologic cycle must on 
the average evaporate as much water from the surface as it precipi- 
tates from the atmosphere in order to conserve water and energy. 
However, the presence of runoff precludes such a balance over land 
alone. That is, the evaporative cooling that is a major determinant of 
summertime land temperatures depends not only on precipitation 
but also on how the surface apportions this input into evapotranspi- 
ration and runoff. This apportioning in turn depends on how net 
radiative energy is divided between latent and sensible heat fluxes. 
Past conventional wisdom has largely ignored potential feedbacks of 
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land on the atmosphere, but with increasing understanding from 
modeling simulations of these exchanges, we can begin to appreci- 
ate what might happen. The most common approach has been to 
study the effect of an arbitrary, though perhaps physically moti- 
vated, change in surface energy or hydrological processes. The study 
of the effect of postulated large increases in surface albedo became 
popular following the suggestion of Chamey et al. (1977) that this 
might represent a positive biophysical feedback. A large number of 
other such studies has been reviewed by Mintz (1984). Of these I 
mention only that of Shukla and Mintz (1982). They considered two 
scenarios, an earth with either perfectly wet or completely dry conti- 
nents Their July simulations showed that the latter would have 
continental surfaces on the order of 10°C warmer than if they were 
wet, and that large reductions in continental rainfall might occur. 
Physical changes of continental land surfaces disrupt global energy 
flows by relatively small amounts compared to increasing trace 
gases (Dickinson, 1986). Therefore, their effects on climate are likely 
to be most pronounced on a regional scale. 

Current studies are attempting to formulate more realistic change 
scenarios, though perhaps still hypothetical. I now consider the 
three kinds of such scenarios. 

The climate effects of tropical deforestation have long been a 
source of speculation. However, now that land surface models are 
beginning to plausibly address the role of vegetation in determining 
the surface fluxes of moisture and energy in GCM climate models, 
there is some hope of establishing at least qualitatively reasonable 
conclusions. The Amazon region contains about half of the world s 
tropical forests, and in the last decade humans have been rapidly 
removing these forests. Thus it is useful to focus the question of cli- 
mate effects of tropical deforestation on this region. Water budget 
studies have established that about half of the precipitation in the 
Amazon is supplied by evapotranspiration from the forest. How might 
removal of the forest reduce the evapotranspiration? Would this 
reduction, in turn, reduce the amount of precipitation? What effects 
might changing surface fluxes and precipitation have on global circu- 
lation patterns by analogy with the more thoroughly studied effects 
of anomalies over the tropical oceans? Several studies in the last few 
years have established that large decreases in evapotranspiration 
would result from forest removal. The most recent studies are also 
indicating comparable decreases in precipitation, possibly even 
greater in amount than the reduction in evapotranspiration. 

Another important role of land surface models is in the explo- 
ration of the contribution of land to year-to-year variations in pre- 
cipitation and temperature anomalies. The continental-scale 
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drought of 1988 and the current long-term drought in California are 
practical examples, dependent to some extent on feedbacks between 
atmosphere and surface. The currently most popular hypothesis is 
that such anomalies are initially related to anomalies in patterns of 
rainfall over the tropical oceans related to the El Nino phenomenon. 
Studies are addressing the role of land interactions in this question, 
but it is too early for definitive results. 

The third scenario currently being studied in which land interac- 
tions have an important role is that of global warming in response to 
increasing greenhouse gases. Model simulations that have produced 
midcontinental drying during the summer have found that the con- 
sequent feedbacks (lack of evapotranspiration) further amplified the 
warming. Interactive clouds have been found to have a yet further 
amplification effect: that is, the surface Is further warmed by more 
solar radiation reaching the surface because of fewer clouds. 


The Process of a Land Model 

The atmospheric components of a GCM (Figure 1) provide the 
surface with fluxes of solar and thermal infrared radiation and pre- 
cipitation in the form of rain or snow, and with near- surface values 
of wind vector, air temperature, and humidity. Water conservation is 
imposed by transferring the water applied to the surface either into 
storage by soil reservoirs or into loss by evapotranspiration or by 
runoff. The total radiation absorbed by the surface is balanced by 
emission of thermal infrared radiation, by the latent heat loss asso- 
ciated with the evapotranspiration or by fluxes of sensible heat, and 
by diffusion of thermal energy into the soil. 

The original Manabe bucket model for the above processes was 
supposed to evaporate at the same rate as a wet surface (zero 
canopy resistance) during well- watered conditions, and to hold a 
maximum water level of typically 0.15 m, which corresponds to the 
available soil water, that is, the water in the rooting zone at some 
average field capacity minus that still present at some average wilt- 
ing point. However, since it did not include the process of diffusion 
of water In soils or canopy resistance, its evaporation rates were 
unrealistic both for bare soil (after a very brief period at the rate of 
atmospheric demand) and for vegetated areas. Evaporation from 
most bare soils, in reality, is greatly reduced after the loss from the 
surface layer of about 1 cm of water. 

Vegetation acts as a completely wet surface only during and 
immediately following precipitation when its foliage is wet. Other- 
wise, it has two important controls: (1) It can extract soil water from 
a greater depth than would evaporation from bare ground, and (2) it 
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atmosphere 



Figure 1. Sketch of the inputs needed at the surface from the 
atmospheric model (Dickinson et al, 1991). 


retards the rate of evapotranspiration from the potential rate 
through resistance by the stomates to molecular diffusion of water. 

The initial bucket model formulation, neglecting any reference to 
canopy resistance, is still widely used in GCMs for studies of climate 
change, and in particular has been used to address the question of 
effects of increasing carbon dioxide and other trace gases (e.g., Man- 
abe and Stouffer, 1980; Washington and Meehl, 1984). The stan- 
dard Goddard Institute for Space Studies (GISS) GCM has a more 
detailed but otherwise similar bucket model with an additional layer 
and water storage capacities depending on ecosystem type, and 
includes a rough approximation of the vegetative removal of deeper 
soil water by allowing infinite upward diffusion during the growing 
season in vegetated regions (Hansen et al., 1983). 
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Inclusion of geographic distributions of soil and vegetation prop- 
erties in a model as shown in Table 1 allows increased realism. In 
particular the aerodynamic resistance of a shortgrass vegetation Is 
over an order of magnitude less than that of forest vegetation, and 
this difference can have major effects on the nature of evapotranspi- 
ration in model simulations. Smooth surfaces have higher tempera- 
tures for the same atmospheric conditions because larger tempera- 
ture and moisture differentials are needed to drive the fluxes 
required by the energy balance (Figure 2). These warmer tempera- 
tures reduce net radiation (by increasing longwave emission) by up 
to several tens of W/m 2 . A comparable additional reduction in net 
radiation is implied by the higher albedo that usually applies to 
shorter vegetation. 

For fixed net radiation and atmospheric conditions, the differ- 
ences in evapotranspiration between tall and short vegetation can 
be obtained from the Penman -Monteith equation. For diy condi- 
tions, the relative effect of changing surface roughness depends on 


Table la: Vegetation /land cover parameters used for South America 





Land Cover/ Vegetation 

Type* 



Parameter 

1 

2 

5 

6 

7 

11 

17 

18 

19 

Maximum fractional 
vegetation cover 

0.85 

0.80 

0.80 

0.90 

0.80 

0.10 

0.80 

0.80 

0.80 

Difference between maximum 
fractional vegetation cover 

0.6 

0.1 

0.3 

0.5 

0.3 

0.1 

0.3 

0.2 

0.3 

and cover at temperature of 
269 K 










Roughness length (m) of 
vegetation 

0.06 

0.02 

0.8 

2.0 

0.1 

0.1 

0.1 

0.8 

0.05 

Depth of the total soil layer (m) 

1.0 

1.0 

2.0 

1.5 

1.0 

1.0 

1.0 

2.0 

1 0 

Depth of upper soil layer (m) 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0 1 

Rooting ratio (upper to 
total soil layers) 

3 

8 

10 

12 

8 

8 

5 

10 

10 

Vegetation albedo for 
wavelengths <0.7 pm 

0.10 

0.10 

0.08 

0.04 

0.08 

0.17 

0.08 

0.06 

0.08 

Vegetation albedo for 
wavelengths >0.7 pm 

0.30 

0.30 

0.28 

0.20 

0.30 

0.34 

0.28 

0.24 

0.30 

Minimum stomatal resistance 
(s/m) 

150 

250 

250 

250 

250 

250 

250 

250 

250 

Maximum LAI 

6 

2 

6 

6 

6 

6 

6 

6 

6 

Minimum LAI 

0.5 

0.5 

1.0 

5.0 

0.5 

0.5 

1.0 

3.0 

0.5 

Stem (and dead matter) area 
index 

0.5 

4.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

Inverse square root of leaf 

10 

5 

5 

5 






dimension (nr 1 / 2 ) 




O 

O 

D 

5 

5 

Light sensitivity factor (W/m 2 ) 

0.01 

0.01 

0.03 

0.03 

0.01 

0.01 

0.01 

0.03 

0.01 


*See definitions in Table lb. 

From Dickinson and Henderson -Sellers, 1988. 
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Table 1 b: Vegetation /land cover types used in the CCM 


1. Crop 

2. Short grass 

3. Evergreen needleleaf tree 

4. Deciduous needleleaf tree 

5. Deciduous broadleaf tree 

6. Evergreen broadleaf tree 

7. Tall grass 

8. Desert 

9. Tundra 

10. Irrigated crop 

11. Semi-desert 

12. Ice-cap/glacier 

13. Bog or marsh 

14. Inland water 

15. Ocean 

16. Evergreen shrub 

17. Deciduous shrub 

18. Mixed woodland 

19. Impoveri shed scrub-grassland* __ 

♦Land type 19 was introduced especially for a deforesta- 
tion study and is not part of the CCM land data set. 


the magnitudes of the net radiation flux and the vapor pressure 
deficit. Under wet conditions the greater surface roughness of 
forests tends to enhance the evaporation loss resulting from inter- 
ception, but how much depends on the spatial distribution of rain- 
fall and on the interception model in use. In the limit of sparse or 
absent vegetation, evaporation will be determined by the treatment 

of diffusion of water through the soil. 

To help summarize the common content of two canopy models 
developed for application in GCMs (i.e., the biosphere-atmosphere 
transfer scheme— BATS— Dickinson et al., 1986, and the simple 
biosphere model-SiB-Sellers et al., 1986), I consider the “lowest 
common denominator” that they both contain, then describe why 
various models introduce further complexity in various features. 
The following derivations (Dickinson et al.. 1991) capture the 
essence of all the model treatments under full canopy conditions 
while leaving out many details such as treatment of fluxes into the 
ground and through leaf boundary layers. 

In meteorological models, the upward flux F x of a quantity X is 
generally represented with the aerodynamic expression 

F x = Pa^"D u (Xs — X a ) ^ 
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Figure 2. Difference between canopy (TJ and atmospheric (TJ temperature 
vs. net radiation (RJfor two different surface resistances (r a ) calculated 
using the simple model described in this paper. Both net radiation and 
surface resistance are defined for isothermal conditions, that is, canopy and 
air temperature the same. Hence radiative and stability feedbacks are not 
included. A relative humidity of 0.7 and air temperature of 22° C are 
assumed; canopy resistance is 100 s/m (Dickinson et at, 1991). 


Subscripts s and a refer to surface and overlying air concentrations 
of X, u is the magnitude of wind, p a the air density, and C D a nondi- 
mensional transfer coefficient. The factor C D is from Monin-Obukhov 
similarity theory for the surface mixed layer of the atmosphere. 

For better semblance to the current notation of micrometeorology , 
we introduce surface resistance r a - Then the flux H of sen- 

sible heat is given by 

H = p a C p {T a -T c )/r a (2) 

where T refers to temperature, C p the specific heat of air, and the 
subscript c refers to the surface being a canopy. For vegetated sur- 
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faces, inclusion of the diffusive resistance by stomates to evapotran- 
spiration E is crucial and so must be included. The integrated effect 
of the resistances of individual leaves is the canopy resistance r c 
such that 

£ = p a (<?a-qc) {3) 

r a +r c 

where q refers to water vapor mixing ratio and q c is determined for 
the internal leaf tissues, i.e., for saturated conditions. The fluxes 
defined by Equations (2) and (3) are illustrated in terms of resis- 
tances in Figure 3. 

Equations (2) and (3) are constrained by the requirement that 
sensible plus latent energy flux be balanced by net radiation R n 
given by 

R n = S^(l - albedo) + Rj± - ecT 4 (4) 


T a’ U 



Figure 3. Schematic resistance diagram for 
the simple generic canopy model described 
here . See text for definitions of terms (Dickin- 
son et al, 1991)., 
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where Sj, = the incident solar energy, Rji = the downward thermal 
infrared energy minus any that is reflected, e = surface emissivity, 
and a = the Stefan-Boltzmann constant. The surface energy balance 
with soil heat flux neglected is written 

H + XE = R n (5) 

where X = latent heat of evaporation. More realistic models add to 
Equation (5) a soil heat flux term. 

Observational data are conveniently analyzed by the Penman- 
Monteith approach of combining Equations (2) and (3) with (5), 
expressing q c in terms of T a - T c and saturated q a (= q a sat ). The 
resulting expression relates E to R n , q a - q c sat > r a , and r c . If all but 
one of these are measured, the remaining quantity, usually either r c 
or E, can be inferred. The main difference between applying these 
principles to analysis of observational data and applying them to cli- 
mate modeling is that for the latter none of these quantities can be 
assumed measured. Rather they must be determined from more 
basic data and model processes. 

Net radiation as described by Equation (4) is conceptually rela- 
tively simple. The atmospheric model provides Sj, and Rji and the 
canopy model provides albedo (according to some combination of 
specified parameters and a canopy radiative model) and T c , the lat- 
ter of which is unknown until Equations (2) through (5) have been 
solved. As already mentioned, r a is obtained from boundary-layer 
theory. The difference q a - q c sat is known from the atmospheric 
model, so the only real complexity is in the specification of r c and 
the determination of T c . 

The leaves are assumed to contribute in parallel so conductances 
l/r s are averaged, i.e., 

r c ~< r s >/LAI (6) 

where the wedge brackets denote an inverse average over the range 
of the canopy leaf area index (LAI), that is, equivalently, the resis- 
tances are summed in parallel. 

We represent the dependence of r s on model variables and for dif- 
ferent ecosystems by a minimum value r smin and a product of limit- 
ing factors, i.e., 

r s - r smln Si (' T)J 2 ( vpd)f 3 {PAR) ... (7) 

where each of the Js has a minimum value of 1 ; /j gives a depen- 
dence on some characteristic temperature, the most obvious being 
that of canopy or root zone soil; f 2 a dependence on vapor pressure 
deficit vpd = (q c - q ca )p s / 0.622 where q ca is the water mixing ratio in 
the air outside the leaves and p s the surface pressure; f 3 a depen- 
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dence on the photosynthetically active radiation flux density (PAR). 
Additional dependences Include water stress (crucial but imple- 
mented differently in different models) and nutrient stress (not yet 
included in existing models). 

Water loss and C0 2 uptake by plants are obviously linked by 
their sharing of the stomates as their path of dominant diffusion 
resistance. Physiologists indicate that this linkage is active; that is, 
stomates act to maintain a constant ratio between water loss and 
carbon assimilation. Such an active linkage between transpiration 
and carbon assimilation would provide a basis for understanding 
the functioned dependences in Equation (7). The calculation of car- 
bon assimilation may be required to determine stomatal resistance, 
or at least may be feasible with little or no additional computation. 
Such a computation would couple the components of the carbon 
cycle with fast time scales to the climate model. 

There is little guidance for current treatments of the temperature 
dependence term, beyond the recognition that optimality will gener- 
ally be achieved in the range of 20-30°C and that stomates will 
cease functioning at temperatures of freezing, 0°C, and of rapid pro- 
tein denaturation, ~50°C. Hence BATS makes a quadratic fit to 
these limits for (Figure 4a). 

There Is also no systematic basis for specifying a vpd dependence 
for r s . However, many observations indicate a near linear depen- 
dence of j >2 1 on vpd , with stomatal closure in the range 0.03 to 0.05 
p s . Hence, SiB and the latent version of BATS assume f 2 l = 1 - 
vpd/c where c 2 0.04 p s (Figure 4b). 

In SiB and BATS, a canopy light model is used to provide light 
levels at a given depth in the canopy, and hence average or integrate 
the/3 component of r s as indicated in Equation (6). This term differs 
superficially between BATS and SiB but is functionally the same 
(Figure 4c). In SiB, it is written 


j - 1 + a 2i c 2 
JJ b 2 + PAR 

where a 2 , b 2 c 2 are adjustable constants. BATS uses 

1 + PAR/PAR C 

r smln / r smax PAR/ PAR c 


( 8 ) 


(9) 


where r smax is the maximum (cuticular) resistance of green leaves, 
and PAR C chosen as the light level where r s = 2 r smin . Equations (8) 
and (9) are equivalent, provided 


r smln ~ c 2 


J It 
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Temperature 



c 





Figure 4. Environmental dependences of inverse of stomatal resistances (le., conduc- 
tance) in BATS model: (a) dependence of conductance on temperature; (b) dependence of 
conductance on vapor pressure deficit, (c) dependence of conductance on photosyntheti- 
cally active radiation flux density, or PAR (Dickinson et al. f 1991). 


— ^2 
r smax ~ c 2 + T 

b 2 

PAR c =b 2 + — 

c 2 

In W/m 2 of visible radiation, PAR,, = 10-50. Precise specification 
of this parameter is neither necessary nor practical. However, fac- 
tor-of-two variations change r c significantly, so ideally, an accuracy 
of better than ±20% is desirable. 

For canopies with random leaf angle distributions exposed to 
direct radiation, the average PAR on a leaf surface is determined from 
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PAR = GrS^ exp(-GD/[i) ( 1 0) 

where D Is the depth Into the canopy in units of LAI, p = cosine of 
solar zenith angle, G ~ 0.5 is the average leaf projection in the direc- 
tion of the sun, and r = ratio of PAR to total incident solar radiation. 
BATS uses Equation (10) with diffuse sky radiation accounted for 
through an additional term which assumes p = 0.5. SiB does like- 
wise but allows in addition for effects of leaf orientation. Neither 
model attempts to account in detail for radiation scattered by the 
foliage, which for PAR has about 0. 1 the intensity of the incident 
solar radiation. SiB uses elegant analytic solutions to Equation (9), 
whereas BATS uses somewhat simpler and probably equally effec- 
tive numerical solutions. 

Interception, the water from precipitation that evaporates from 
the canopy without reaching the soil, has similar one-layer parame- 
terizations in all models. The interception parameterizations also 
provide dew or frost formation when the water vapor gradient from 
foliage to air reverses. The fraction of canopy surface covered by 
water has zero resistance r c and hence, especially In forests, can 
rapidly evaporate back into the atmosphere. The parameterization 
choices are the water-holding capacity of the foliage W sc , the frac- 
tion of the incident precipitation that Is intercepted, and the fraction 
of the foliage that is covered by water when it has less water than 
full capacity. After the canopy reaches capacity, all additional pre- 
cipitation is put into throughfall. 

BATS and SiB now use W sc = 0.1 LAI (in mm) for water capacity. 
This quantity is tuned to give observed canopy interception losses 
and is somewhat lower than the observed storage of water by leaves. 
BATS has a “stem area index” surface, as well as LAI, that is also 
wetted. SiB determines a cross section for interception depending on 
LAI, similar to that for radiation, whereas BATS assumes all precipi- 
tation over vegetation is first captured by foliage. SiB assumes the 
fraction of canopy surface wetted is the ratio of canopy water to W sc , 
whereas BATS and a new model being developed at GISS use a 2/3 
power. In reality, the fractional wetting is veiy dependent on the 
hydrophilic properties of the leaves. Some leaves are partially wetted 
with smooth water films, while on others, water droplets form. Thus, 
for application to specific sites and vegetation, the present models 
may be unrealistic. The values of W sc in SiB and BATS are inferred 
from two years of Amazon measurements (Shuttleworth, 1988). 

Model-Specific Parameterizations 

The components of vegetation resistance summarized in the pre- 
vious section either are treated similarly in essentially all of the veg- 
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etation resistance models or else are assumed constant. While 
details are debatable, we can agree that there is an appropriate 
functional form consistent with observational information, and that 
the prescribed constants in the different models are not drastically 
different. In this section I discuss those components that diverge 
more drastically in the different land parameterizations. 

Water Stress 

The models differ in their treatment of the effect of water stress on 
stomata! resistance, In part because of the lack of reliable quantita- 
tive information on the subject. SiB assumes that this contribution 
to r s depends on leaf water potential, with the leaf water potential 
being related to soil water potential through the effect of soil and 
root diffusion resistance to the water movement. These processes 
have been represented through more detailed mechanistic models, 
and some observational information is available for individual sites. 
However, there is little or no basis for specifying the necessary para- 
meters over large areas. 

BATS uses a simplified version of the approach illustrated in Fig- 
ure 5. Each ecosystem is characterized by a maximum transpiration 
rate under well-watered conditions, presumed to be determined by 
root and soil resistance. This maximum rate is reduced below field 
capacity according to the difference between the soil water suction 
(negative potential) for wilting and that computed for existing soil 
moisture. The latter is obtained as an average over the model soil 
layers, where the average is weighted with the root surface density 
in each layer. This averaging requires only an estimate of the rela- 
tive distribution of roots because their absolute contribution is sub- 
sumed in the assumed maximum transpiration rate. The most cru- 
cial parameter in the BATS treatment of root resistance is this 
maximum rate, which in principle could be specified from remote 
sensing. 

A similar but even simpler concept is used in the new GISS model 
being developed, i.e., the contribution to transpiration is reduced in 
each soil layer from that computed by the unstressed canopy model 
by the ratio of the total water potential to that for wilting. 

Within-Canopy Resistances 

SiB and BATS attempt to include the bulk effect of boundary- 
layer resistances across leaf surfaces. Given a local wind in the 
canopy, this resistance for heat and moisture can be inferred from 
laboratory studies. Its accurate specification is limited by the knowl- 
edge of the wind distribution within the canopy. BATS simply esti- 
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Figure 5 , Schematic 0 / the approach used by BATS to 
determine dependence on stomatal resistance on soil 
water. Atmospheric demand is determined and com- 
pared with the maximum that can be supplied by roots 
for given soil water. If the former exceeds the later , 
stomatal resistance is increased to reduce the demand 
to match supply (Dickinson et al, 1991). 


mates this wind from the frictional velocity, whereas SiB bases it on 
a solution for eddy diffusion within the canopy. 

Partial Vegetation 

The treatment of partial vegetation in BATS helps illustrate some 
of the questions that must be faced in treating this issue. The frac- 
tion of vegetation covered by a model grid square is prescribed, with 
a seasonal variation determined from soil temperature. No specifica- 
tion is given as to the spatial scales to be associated with the bare 
soil fraction. Land classes range from desert (that is, all bare soil) 
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and semidesert (that is, mostly bare soil) to various forest types that 
are mostly vegetation cover. Sensible and latent fluxes are computed 
separately for the bare ground and for the fraction of soil under veg- 
etation. Both soil fractions are assumed to have the same tempera- 
ture and moisture content but to differ in values of overlying wind 
and transfer coefficients, i.e., different values of r a are determined 
for these two fractions. The details of these prescriptions are guided 
by their reasonableness in the limits of bare soil and full canopy. A 
more realistic treatment might determine separate soil temperatures 
and moisture for shade vs. sun and require information on the spa- 
tial scales of the vegetation and exposed soil areas. 

In SiB, the effects of partial vegetation are incorporated directly 
into the radiation, momentum (turbulent), and energy transfer sub- 
model. However, no account is taken of larger- scale heterogeneities. 
The new GISS model separates the grid box into vegetation-covered 
and bare soil components. Fluxes and soil temperatures are calcu- 
lated for each of these surface types, and then area-weighted to 
interface to the first layer of the atmosphere. 

Canopy Temperature 

Calculation of canopy temperature has been the most difficult 
aspect of vegetation resistance models in GCMs, because implement- 
ing the calculation of surface temperatures can lead to worrisome 
inaccuracies and, at worst, to severe computational instabilities. 
Without a successful approach to this question, much greater errors 
in determination of evapotranspiration can be made than might 
result from inaccuracies in canopy resistance. In general, we must 
distinguish between two or more temperatures that are coupled 
through energy fluxes and must separately satisfy energy balance 
requirements. For example, SiB and BATS have separate soil and 
canopy temperatures, and the GISS model in addition distinguishes 
between bare ground and under- vegetation soil temperatures. These 
surfaces each have heat capacities, some of which may be zero. 

Let T be the vector representing all the model surface tempera- 
tures and C a diagonal matrix where elements are individual surface 
heat capacities. To determine T, knowing its value at a previous time 
step, we must solve numerically an equation of the form 

c f ~ F(T)=0 ai) 

where F(T) is a vector whose individual elements represent the sum 
of energy fluxes into a given surface and t is time. Let superscript n 
refer to the value of T at the nth model time step. Time steps are At 
We could first try the simplest solution to Equation (11), i.e., 
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T n =T rl - 1 +C- 1 F(T n " 1 )A£ (12) 

This may work with short enough time steps, a few minutes or less, 
and large enough heat capacities (i.e., provided At is small compared 
to the inverse of the largest eigenvalue of C -1 5F/5T), but otherwise 
it can be a prescription for numerical disaster giving wild and grow- 
ing oscillations from one time step to the next. In the limit of small 
heat capacities, Equation (11) should approach a statement of 
energy balance at the present time level. An alternative solution 
likely to be more accurate and stable is hence 

T a =T n_1 +C _1 F(T ri )At (13) 

but this form appears to require already knowing the solution to 
use. However, its solution may be possible using a Newton iteration, 
i.e., by writing the fth component of F = F j as 

Fj (T".^ 1) = Fj(Tn-t)+ SFj jj n ’ --(T^’ <+1 -TQ*) (14) 

0 1 k 

where the second, ith, superscript refers to the number of the itera- 
tion, there is a summation over the kth subscripts, and T k is the kth 
temperature. Equation (14) is substituted into Equation (13) for 
FfP V and Equation (13) is solved for T n t+1 , taking as a first guess 
the value of T at the previous time step, that is, T n l = T" -1 . 

The SiB model specifies soil and canopy heat capacities and uses 
only the first-guess form of Equation (14). BATS, on the other hand, 
assumes zero heat capacity for the canopy and first iterates Equa- 
tion (14) for canopy temperature to convergence, taking soil temper- 
ature as that from the previous time step. It then solves Equation 
(13) as a scalar equation for soil temperature, using the first derived 
canopy temperature and only the first guess from Equation (14). A 
further simplification, sometimes used, e.g., at the European Centre 
for Medium Range Weather Forecasts (ECMWF), is to assume all 
surfaces have the same temperature, so that Equation (13) is 
applied as a scalar to derive joint soil and canopy temperatures. 

Preliminary examination of the errors from the SiB and BATS 
solutions for canopy temperature, including lack of conservation of 
energy, suggests that both approaches can reasonably control 
errors; SiB conserves energy even if its first iteration is inaccurate, 
whereas BATS only does so with convergence of the temperature 
iteration. The BATS solution may be more accurate for canopy-dom- 
inated transpiration, but the SiB approach requires less numerical 
computation and so may be preferable, considering all of the other 
uncertainties in the parameterizations. 
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Present Experience with a 
L an d Model Coupled to a GCM 

Here we will report on the latest simulation studies with the com- 
munity climate model version 1 (CCM1) of the National Center for 
Atmospheric Research (NCAR) coupled to the BATS land surface 
treatment. Shortcomings are emphasized more than successes to 
point to where further progress is needed. 

Over the last year, we have developed two slightly different versions 
of the CCM 1 tuned to give top-of-the-atmosphere fluxes of solar and 
longwave radiation that match those observed by satellite. Both start 
with the standard CCM1 with a diurnal cycle added and the solar 
radiation calculation done with a treatment developed by A. Slingo of 
NCAR and the British Meteorological Office. One version in addition 
calculates clouds with a scheme developed by A. Slingo, initially for 
the ECMWF model. During these studies, we discovered that previous 
models which assumed a diurnal average sun would absorb globally 
about 10 W/m 2 additional solar radiation when a diurnal cycle was 
added, hence requiring us to adjust the solar cloud treatment to 
reflect about that much more radiation to get the match with data 
from the Earth Radiation Budget Experiment (ERBE) satellite (since 
the treatment had previously been timed for an average sum). 

The bulk of our simulations have been done with the standard 
CCM1 clouds, and we have not noticed any major differences over land 
with the other code. The integration includes an eight-year control sim- 
ulation with a lower boundary over ocean as the new sea surface cli- 
matology developed by K. Trenberth and collaborators. The first four 
years of this integration are used to specify monthly average ocean 
heat transport (flux correction) needed to balance net surface energy 
fluxes at each ocean grid point. Otherwise, we treat the ocean at each 
grid point as a 60-m slab of water for its temperature calculation. 

Together with a thermodynamic sea ice model, we have begun to 
integrate this CCMl/BATS/flux-corrected slab-ocean model for 330 
and 660 ppm of C0 2 . In both cases we used a 30-m slab for one 
year to speed up convergence and, for the doubled C0 2 , initialized 
with a uniform 3°C warming over all the ocean points. 

The purpose of these integrations is to explore the behavior and 
feedbacks of a detailed land model (BATS) in a global climate change 
scenario, comparable to those that have been used previously with 
simple bucket models. 
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